Bose-Fermi solid and its quantum melting in an one-dimensional optical lattice 
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We investigate the quantum phase diagram of Bose-Fermi mixtures of ultracold dipolar particles 
trapped in one-dimensional optical lattices in the thermodynamic limit. With the presence of 
nearest-neighbor (N.N.) interactions, a long-ranged ordered crystalline phase (Bose-Fermi solid) is 
found stabilized between a Mott insulator of bosons and a band-insulator of fermions in the limit 
of weak inter-site tunneling (J). When J is increased, such a Bose-Fermi solid can be quantum 
melted into a Bose-Fermi liquid through either a two-stage or a three-stage transition, depending 
on whether the crystalline order is dominated by the N.N. interaction between fermions or bosons. 
These properties can be understood as quantum competition between a pseudo-spin frustration 
and a pseudo-spin-charge separation, qualitatively different from the classical picture of solid-liquid 
phase transition. 

PACS numbers: 



The successful experimental preparation of ultracold 
dipolar atoms []J and polar molecules [2| has opened a 
new direction of strongly correlated physics. The long- 
ranged anisotropic feature of the dipolar interaction can 
lead to qualitatively different physics as compared to the 
conventional ultracold atom systems. Recent theoreti- 
cal efforts have been mostly devoted to either bosonic or 
fermionic dipolar gases [3 5]. However, not much work is 
reported on dipolar Bose-Fermi mixtures, although their 
counterparts in conventional ultracold atoms have been 
extensively studied both experimentally [6| and theoret- 
ically [7H10j. When the dipolar interaction is included, 
some crystalline order may be realized Q , providing com- 
pletely new insight onto the solid-liquid phase transition, 
which is one of the most important subjects in classical 
physics. We show in this work that a Bose-Fermi dipo- 
lar mixture allows for extremely rich interaction physics 
by manifesting a complex and experimentally accessible 
quantum phase diagram containing many exotic quan- 
tum phases hitherto not considered in cold atomic gases. 

In this paper, we study various kinds of Bose-Fermi 
mixtures of dipolar particles loaded in an one dimensional 
(ID) optical lattice. Keeping only the on-site and the 
nearest-neighbor (N.N.) interactions, five distinct quan- 
tum phases can be identified: Bose-Fermi solid (BFS), 
density wave Bose-Fermi Mott state (DW-BFM), den- 
sity wave Bose-Fermi liquid (DW-BFL), uniform Bose- 
Fermi Mott state (UBFM), and uniform Bose-Fermi liq- 
uid (UBFL), as shown in Fig. 1. The BFS phase has al- 
ternating boson and fermion densities similar to its clas- 
sical counterpart, but can be quantum melted even at 
J = through an Ising-type frustration when the dipole 
moments of bosons and fermions are equal. When J in- 
creases, the BFS state can have another quantum melting 
through a multi-stage transition, which can be under- 
stood as Kosterlist-Thouless type with the pseudo-spin- 
charge separationjllj. We also discuss the experimental 
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FIG. 1: Phase diagrams of dipolar Bose-Fermi mixtures: (a)- 
(d) corresponds to Case I-IV as defined in the text. Vbb —0.1 
in (b) and (d), and Vff = 0.225 in (c) and (d). Vbf = 0.15 
in (d) and /if = 0.48 for all panels. U = 1 is the energy 
unit. The vertical dotted lines in (a)-(c) show the intersect 
for which the particle density is presented in Fig. 3. 



preparation and observation of these quantum phases, 
and connect our results to the metallic electronic recon- 
struction at the interface of a band insulator and a Mott 
insulator [l2| . 

To describe our model system, we adopt the following 
extended Bose-Fermi Hubbard Hamiltonian: H = Hb + 
Hp + Hbf, where 
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Here, bi and fa are the bosonic and fermionic field opera- 
tors at the zth site, nf = b\bi and nf = fjfi are density 
operators. J a and fi a are respectively the hopping rates 
between adjacent sites and chemical potentials; 

U a b and V a b describe the on-site and the N.N. interac- 
tions respectively, with a and b denoting bosons (B) or 
fermions (F). Such a system can in principle be engi- 
neered in ultracold dipolar gases. Longer range interac- 
tions are neglected here because they do not qualitatively 
change the phase diagram near half filling 14j. 

In the hard-core boson limit (Ubb — > oo) one can 
map the bosonic operators to fermionic ones, and the 
whole system can be transformed to a special extended 
Hubbard model with spin-dependent N.N. interactions. 
Such a system breaks the SU(2) symmetry and cannot 
be solved analytically. Since we are interested here in 
the interplay of quantum statistics with N.N. interac- 
tions, we will consider a soft-core case by setting J a = J 
and U a b = U = 1 for simplicity. Near the Mott insu- 
lator regime, we can reduce the Hilbert space into 0, 1 
and 2 particles per site and hence such soft-core bosons 
can also be mapped to a spin-1 chain [U, i.e. Eq. (p} 
can be mapped onto a unique spin ladder composed of 
a spin-one chain and a spin-half chain with frustrated 
anti- ferromagnetic coupling as V a b > 0. In our calcula- 
tion, we stick with the particle model and the ground 
states are obtained using an infinite lattice version of the 
time-evolving-block-decimation (TEBD) algorithm [l3j ]. 
which exploits a tensor product representation of many- 
body states and is especially efficient for the gapped solid 
phase that we are interested in here. 

In order to clarify the effects of the N.N. interactions, 
we calculate and compare the following six representa- 
tive cases: Case I: neither bosons or fermions are dipolar 
particles, i.e. V a b = 0; Case II: only bosons are dipolar 
particles, i.e. Vbb > and Vbf = Vff = 0; Case III: 
only fermions are dipolar particles, i.e. Vff > and 
Vbb = Vbf = 0; Case IV: both bosons and fermions 
are dipolar particles, i.e. Vbb /Vbf = Vbf /Vff = X 
and x = 2/3 < 1; Case V: same as Case IV but 
for x = 3/2 > 1; Case VI: same as Case IV but for 
X = 1- Several representative quantum phase diagrams 
are shown in Fig. 1, where panels (a)-(d) correspond to 
Case I-IV respectively. Results of Case V and Case VI are 
qualitatively the same as Case II and Case I respectively. 
To save space, we do not show their phase diagrams here. 

In the limit of small (large) /xb, the system is a pure 
Fermi(Bose) gas, denoted by F(B). Within these param- 
eters, the ground state can be a bosonic Mott insulator 
or a ID Bose liquid in regime B, while in regime F it 
can be a band insulator or a Fermi liquid. In the mixed 
regime, there are five distinct quantum phases to be ob- 
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FIG. 2: Typical real space distributions of boson (open trian- 
gle), fermion (filled circle), and total density (open diamond) 
in (a) UBFL, (b) UBFM, (c) BFS, (d) DW-BFM, and (e) 
DW-BFL phases. Panel (f) shows a BFS phase with higher 
filling fraction of bosons. 



served: (i) in the large J regime, we have a UBFL phase, 
whose quasi-long-ranged order can be obtained by us- 
ing bosonization[8]. The density distribution is uniform 
in space for each species, as shown in Fig. 2(a). (ii) 
In the regime of stronger interaction, there is a UBFM 
phase, where the total filling fraction is found to be unity 
through the entire regime with a gapped excitation. The 
UBFM phase exists even without N.N. interaction (Fig. 
1(a)) [101, and it has a uniform and compressible inter- 
species density difference (Fig. 2(b)). (iii) When the 
N.N. tunneling is close to zero, we find a BFS phase, 
where both species are at half-filling but have an alter- 
nating density pattern (see Fig. 2(c)). Such a phase is 
incompressible for both species (i.e. d(nf)/dfib = for 
a, b = {B,F}), while their density oscillation amplitudes 
decrease as J increases, (iv) Near the boundary between 
BFS and UBFM, there is a narrow regime, where the 
BFM phase has a density wave order (DW-BFM), i.e. 
(nf — nf) oscillates periodically while the total filling 
remains integer (see Fig. 2(d)). (v) Finally, near the 
boundary between BFS and BFL, we find a BFL with 
a density wave order in both species (named DW-BFL). 
Different from DW-BFM, DW-BFL is compressible with 
an incommensurate total filling fraction (Fig. 2(e)). This 
phase can also be regarded as a supersolid phase in the 
quasi-long-ranged order sense. 

In Fig. 2(a)-(e), we show typical spatial density dis- 
tributions for BFL, UBFM, BFS, DW-BFM, and DW- 
BFL phases respectively. We note that, although the 
true spin-charge separation exists only in the low en- 
ergy regime [8|, we can still phenomenologically define 
the pseudo-charge to be (nf + nf), and pseudo-spin to 
be (nf ~ nf). As a result, the BFL phase has gapless 
excitations in both pseudo-charge and pseudo-spin sec- 
tors, UBFM has a gapped charge excitation, and BFS 
has gapped excitation in both sectors. DW-BFM has a 
gapped pseudo-charge excitation with a gapless pseudo- 
spin-density wave order, while DW-BFL has two gapless 
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FIG. 3: (a)-(c) Average boson (circle), fermion (triangle), and 
total (diamond) densities as functions of [ib at J — 0.1. All 
other parameters are the same as in Fig. l(a)-(c) respectively, 
(d) shows the density wave order at J = 0.17 for the same 
system as in (c). The solid dot and check markers correspond 
to the density wave order, \{nf) — (nf +1 }\ and |(nf ) — {nf +1 }\. 



excitations on top of the pseudo-spin density wave or- 
der. Note that our numerical calculation is for an in- 
finite ID system, so the observed density wave order is 
always a truly long-ranged order. In Fig. 3(a)-(c), we 
further plot the average filling fraction of bosons and 
fermions (together with the total fraction) as a function 
of the bosonic chemical potential \ib > showing the regime 
of different phases. The plateau for particle density at 
0.5 in panels (b) and (c) corresponds to the half-filled 
BFS phase. The plateau of total density at unity can 
be identified to be a Mott phase if densities of individ- 
ual species vary as a function of the chemical potential. 
Finally, the regime with a varying total density corre- 
sponds to the BFL phase. UBFM and DW-BFM can- 
not be distinguished here (same as UBFL and DW-BFL) 
since only the average filling fraction is shown. In Fig. 
3(d), the density wave order is plotted as a function of 
Hb at J = 0.17C/, and is found to disappear continuously 
when changing fis away from BFS regime. 

Note that the origin of our BFS state is different from 
that of the Neel(Ising) state predicted in where no 
N.N. interaction exists. In the latter case, the effec- 
tive long-ranged interaction between heavier species is in- 
duced by the quantum fluctuations of the lighter species, 
and becomes effectively vanishing when Jg = Jp. Be- 
sides, we also note that at higher fiB (outside the regime 
of Fig. 1), a BFS state with more than one particle per 
site is also possible (see Fig. 2(f)), but we will not discuss 
such phases here due to space limitation. 

To gain more insight on the crystalline phase, it is in- 
structive to look into the classical limit at J = 0, where 
no difference between fermionic and bosonic particles is 
expected. By comparing the energies of various con- 
figurations, for example, half-filled Bose/Fermi density 



wave, Bose Mott Insulator, and half-filled BFS, we can 
obtain the necessary condition for a BFS phase to be: 
UBiHf > 2Vbf and 

2{V FF - V BF ) > fJ-F - m > 1{V BF - Vbb). (2) 

It is therefore interesting to find that for a symmetric 
case, i.e. Case VI and V a b = V, a crystal exists only when 
= > 2F, and the state {• • • BFBF ■ ■ ■} would be 
degenerate with all states of one particle per site. This 
can also be understood in the spin-ladder picture we de- 
scribed earlier: the N.N. interactions, V a b > 0, provide an 
Ising-type frustration, namely, Vbf competes with Vbb 
and Vff, and hence smear out the anti- ferromagnetic or- 
der. At a finite but small J, the degeneracy is removed, 
but the density wave order is still so small that the BFS 
phase just disappears when [1b slightly differs from /ip- 
Therefore, according to our numerical calculation and 
conditions in Eq. (2), we conclude that the BFS phase is 
unlikely to be observed in a realistic experimental situa- 
tion for the symmetric case (x = 1) due to the frustrated 
nature of the dipolar mixture. The system is therefore 
still compressible and uniform in density, as the UBFM 
case shown in Fig. 1(a). 

Away from the symmetric case, we first consider the 
pro-boson case, i.e. the N.N. interaction between bosons 
is stronger than that between fermions (x > 1) and hence 
dominates the crystal order formation. For such a situa- 
tion, we show our results for Case II in Fig. 1(b). Results 
for Case V are similar. Such pro-boson BFS is quan- 
tum melted to the UBFL via a three-stage process: first, 
the ground state has a second order transition toward 
the DW-BFM phase, having a compressible pseudo-spin 
density wave order on top of a pseudo-charge Mott phase. 
When J becomes slightly larger, the density wave order 
disappears and the ground state becomes a UBFM phase, 
where the single particle Green's function decays expo- 
nentially (see Fig. 4). Finally, the pseudo-charge gap 
of UBFM becomes zero at larger J, leading to a BFL 
phase. The last two transitions should belong to the 
Kosterlitz-Thouless universality class in the pseudo-spin 
and pseudo-charge sectors respectively, as described by 
an X-Y model in the 1+1 dimensional system. 

However, such kind of three-stage process from BFS 
to BFL does not always exist for a pro-fermion case, i.e. 
when the crystal order is mainly provided by the N.N. 
interaction between fermions. We provide the quantum 
phase diagrams for both Case III and Case IV in Figs. 
1(c) and (d). In Fig. 1(c), we find that in addition to 
the three-stage melting process discussed above, in the 
lower portion of the BFS boundary, there is a different 
melting process: first from BFS to DW-BFL and then 
from DW-BFL to UBFL. 

Such different melting processes for pro-boson and pro- 
fermion solids demonstrate the quantum character of 
BFS. Their common feature is the smooth decrease of the 
density wave order when outside the BFS phase, while 
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FIG. 4: Green functions for (a) bosons and (b) fermions in a 
few samples inside the UBFM and UBFL phase. J = 0.15, 
0.16, 0.17, and 0.18 for the solid, dashed, dotted, and dash- 
dotted lines respectively. fj,B = 0.4 and all other parameters 
are the same as in Fig. 1(c). 



the difference fies in that the pseudo-charge gap disap- 
pears together with the pseudo-spin gap in the lower por- 
tion of BFS boundary for the pro-fermion case. This re- 
sult can be understood when comparing the BFS regime 
in Fig. 1(b) and (c). First, without N.N. interaction (Fig. 
1(a)), the BFS is degenerate with other states at a single 
point, hb = an d J = 0. When Vbb > and other 
V a b = 0, from the condition shown in Eq. ©, the regime 
of BFS is < fin - fJ-F < 2Vbb, i-e. the BFS regime is 
entirely inside the original Mott insulator regime. As a 
result, when BFS is melted, the whole system is still in 
the regime of commensurate total filling fraction, and a 
Mott phase with a pseudo-charge gap is then expected 
in the strong interaction limit. On the other hand, when 
only Vff > and other V a b — 0, the condition for BFS 
at J = gives < fiF ~ A's < 2Vff, which is inside 
the regime of fermion band insulator as shown in Fig. 
1(a). Unlike Mott insulator for bosonic particles, the 
band insulator of fermions exists in non-interacting limit 
and hence can be easily destroyed by finite interactions, 
leading to a Luttinger liquid in ID with incommensurate 
filling. That is why both pseudo-spin and pseudo-charge 
gaps disappear in the DW-BFL and UBFL regimes when 
BFS is melted in the lower boundary of BFS phase in 
Fig. 1(c). The nature of such a transition is still unclear. 
Finally, from Fig. 1(d), we can see that the phase dia- 
gram for such a pro-fermion case is qualitatively similar 
to Fig. 1(c), but the finite value of Vbf can suppress the 
BFS phase as well as the BFM phase, consistent with 
the fact that no BFS exists in the symmetric case due to 
Ising-type frustration. 

In Fig. 4 we also present the single particle Green's 
function of (a) bosons: Gfl(r) = [b\bi +r ) and (b) the 
fermions, Gf(t) — (fjfi+r)- It clearly shows that the 
single particle correlation of bosons and fermions has an 
exponential decay in the UBFM regime (solid and dashed 
lines), while they have a slower algebraic decay due to the 
quasi-long-ranged order inside the UBFL regime (dotted 
and dashed-dotted lines). 

Before concluding, we note that the system we dis- 
cuss here can be easily realized by a mixture of dipolar 



atoms/molecules loaded into a ID optical lattice. One 
may provide two separate trapping potentials for bosons 
and fermions, and prepare them to be a Mott insula- 
tor and a band insulator respectively. When the relative 
distance between these two traps is getting smaller, the 
inter-mediate regime between these two clouds experi- 
ences a continuous change of chemical potentials, making 
a potential realization of BFS and other exotic phases dis- 
cussed here. Our results indicate that observing a BFS 
could be very challenging for systems in which bosons 
and fermions have close dipole moments (say 39 K 87 Rb- 
40 K 87 Rb or 52 Cr- 53 Cr), due to strong frustration. It is 
therefore more promising to consider a mixture of dipolar 
and non-dipolar particles to realize a BFS. It is interest- 
ing to note that, compared to the metallic and ferro- 
magnetic interface of an electronic band insulator and a 
Mott insulator [l2[, the intermediate BFS phase here is 
an insulator with anti-ferromagnetic order, which can be 
easily observed in a time-of-flight measurement. 

In summary, we have studied the ground-state phase 
diagram of the Bose- Fermi Hubbard model with arbitrary 
N.N. interactions. We discover the possibility and the 
condition to find a novel Bose-Fermi solid phase, which 
has a nontrivial quantum melting due to the mixture of 
quantum statistics even in ID. This system can be ex- 
perimentally realized within the present technique. 
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